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Abstract 

Recent analysis of social communications among humans has revealed that the interval between 
interactions for a pair of individuals and for an individual often follows a long-tail distribution. We 
investigate the effect of such a non-Poissonian nature of human behavior on dynamics of opinion 
formation. We use a variant of the voter model and numerically compare the time to consensus of 
all the voters with different distributions of interevent intervals and different networks. Compared 
with the exponential distribution of interevent intervals (i.e., the standard voter model), the power- 
law distribution of interevent intervals slows down consensus on the ring. This is because of the 
memory effect; in the power-law case, the expected time until the next update event on a link is 
large if the link has not had an update event for a long time. On the complete graph, the consensus 
time in the power-law case is close to that in the exponential case. Regular graphs bridge these 
two results such that the slowing down of the consensus in the power-law case as compared to the 
exponential case is less pronounced as the degree increases. 

PACS numbers: 02.50.-r, 05.40.-a, 89.65.Ef, 89.75.Fb 



I. INTRODUCTION 



Macroscopic social dynamics often occur as a result of microscopic dynamics of individuals 
interacting on networks of social contacts. Studies of interacting particle systems such as 
spin systems have enriched our understanding of various types of social dynamics such 
as epidemics, information cascades, opinion formation, synchronization, and evolutionary 
games jl-4]. It is established that the structure of social networks influences social dynamics 
in many different ways. 

Social dynamics on networks are often modeled by stochastic processes. Models of these 
dynamics usually assume that interaction events between a pair of individuals occur accord- 
ing to the Poisson process. This assumption facilitates theoretical analysis of models and 
corresponds to the situation where the rate at which an event occurs generally depends on 
the current configuration of the network but not on the history of the dynamics. 

Recent developments of sensor technologies and accumulation of massive amounts of 
electronic data have facilitated detailed analysis of point processes related to social d yna mics 



and face-to-face conversations 



of humans. Examples of such data include email exchanges js 9 ] , cell-phone calls jlOl . Illfl . 



12h14|. Apart from the network structure, these studies 



have revealed a novel universal feature of social dynamics: the non-Poissonian interevent 
intervals (IEIs). The distributions of IEIs are often inherited with long tails and can be 
modeled by the power-law distribution possibly with an exponential cutoff 0, Q] or by 



the log-normal distribution 



prioritization of tasks 
of human activities 



161 ] . The long-tail IEI distribution can be explained by the 



7| or by the combination of seasonality and the circadian rhythm 



The effect of the long-tail IEI distributions on epidemic dynamics such as the suscepti 



)lc- 



infectedfSI) and susceptible-infected-recovered (SIR) models has been investigated 



r] 



18 



3,H 



-l20j. It has been suggested that the long-tail IEI distribution is responsible for the 
persistent prevalence of computer viruses [ijj] and email advertisements 16 1 that cannot 
be explained by the Poisson assumption. The slowing down of epidemic dynamics owing 
to the long-tail IEI distributions is also found in the epidemic model on priority-queue 
networks |19| . 

In this paper, we consider the effect of the long-tail IEI distribution on opinion dynamics. 
We use a variant of the voter model and compare the time required by this variant to 



achieve consensus among all the voters with that required by the standard voter model with 
the exponential IEI distribution. By numerical simulations, we show that the long-tail IEI 
distribution increases the consensus time on the ring. The consensus time obtained for the 
exponential and long-tail IEI distributions is relatively close on the complete graph. We 
interpolate the results on the ring and the complete graph by examining the voter models 
on the regular random graph; by this interpolation, we show that the node degree is a main 
determinant of the difference in the consensus time. 



II. MODEL 



We analyze a variant of the voter model [2|, l2lH23| on static regular networks; i.e., all 
the nodes have the same degree. A voter is placed at each node in the network, and each 
voter takes one of the two opinions denoted by and 1. Initially, each voter takes 
and 1 with equal probability (i.e., 0.5), and the opinions of different voters are assigned 
independently. Each link between the nodes is independently endowed with a random IEI 
Ti, which represents the time until the initial update event occurs on this link. We denote 
the distribution of T\ by Pi(ti). Suppose that an update event on a link occurs at a certain 
time. If the two endpoints of the link are occupied by the opposite opinions, one of the 
two voters is selected with equal probability (i.e., 0.5) and the opinion of the selected voter 
is flipped such that the two voters take the same opinion. Otherwise, nothing happens in 
the update event. Then, the next IEI for this link is drawn from distribution p(r). The 



sequence of update events on each link is a renewal process [24( and only the initial IEI T\ 
obeys p\(r\), whereas all the subsequent IEIs obey p(r). If T\ were taken from distribution 
p(r), we would be implicitly assuming that there is an event at time 0. Introduction of Pi(ti) 
is necessary to avoid such artificial initial conditions. Update events and the assignment of 
random IEIs occur independently on all the links. It should be noted that update events 
never (i.e., with measure zero) occur on multiple links at the same time if Pi(ti) and p(r) 
do not possess point mass, which we assume. We assume that each node is updated once 
per unit time such that 6t = 1/N. 

In this study, we assume that Pi(ti) obeys 

p 1 (T 1 ) = ^- ) J~p(T>)dT', (1) 



where (•) denotes the mean. With the definition of P\{t\), the renewal process is the so-called 



equilibrium renewal process 



24|. If the renewal process is ongoing and the voter dynamics 



begin at an arbitrary instant, then T\ obeys P\[j\) given by Eq. ([I]), not p(r). 

We run the dynamics until the entire network is taken over by one opinion; we refer to 
such a unanimous configuration as a consensus. We are concerned with the consensus time, 
i.e., the time required to reach a consensus. To obtain the averaged quantities, we perform 
1000 rounds of simulations under each condition unless otherwise stated. 

In the present study, we examine the following two voter models. 



For the exponential voter model, we set 



p(r) = Aexp(-Ar), 
Pi(ti) = Aexp(-Ari). 



(2) 
(3) 



For the power-law voter model, we set 

a — 1 / r + c 



p(r) 



a - 2 



Ti + C 



(4) 
(5) 



In the exponential voter model, the occurrence of events on each link obeys the independent 
and identical Poisson process. The exponential voter model is equivalent to the standard 
voter model with the updating procedure called link update [25] or link dynamics (LD) 



26|. 



For the power-law voter model, our choice of the power-law p(r) is motivated by recent 



experimental results obtained for face-to-face interactions 



the power-law p(j) has a finite mean. Empirically, p(r) usua 
two |l3-8, 10-12, 15 1 and accompanies an exponential cutoff 



121 ] . We assume a > 2 such that 
ly has a value of a less than 



7 



10 



15j. With a < 2, the 



consensus time of the power-law voter model trivially diverges for any network because (r) 
diverges. Instead of setting a < 2 and modulate a or the cutoff value of r, we use Eq. (Tj|]) 
with different values of a > 2 to examine the effect of the large variance of interevent 
intervals on the consensus time. 

We set A = 1 and c = a— 2 to set (r) to unity for both the exponential and power- law voter 
models. In our numerical simulations, we set (a, c) = (2.5,0.5) and (3.5,1.5). We verified 
that our numerical results presented in the following sections for a = 2.5 are qualitatively 



the same in the range 2 < a < 3 and that the results for a = 3.5 are qualitatively the same 
in the range 3 < a < 4. 

The heterogeneous degree distribution of the network, which is eminent in scale-free 



networks, makes the consensus time sensitive to the adopted update rule [25|, [26|. The 
consensus time in networks with heterogeneous degree distribution has been investigated for 
the exponential (i.e., standard) voter model 25j-|3l|. To focus on the effect of the power-law 
IEI distribution, we restrict ourselves to the regular graphs in the present study. Specifically, 
we compare the consensus time of the exponential and power-law voter models on the ring, 
the complete graph, the extended ring, and the regular random graph. 



III. RESULTS 



Ring 



Assume that the voters are placed on the ring with N nodes. The degree of each node is 
equal to two. 

To understand the mechanism governing the consensus time, denoted by T, we begin with 
tracking the fraction of voters who take opinion 1 and the number of the interfaces, denoted 
by m and E^, respectively. For the ring, a link is defined to be an interface when the two 
endpoints of the link are occupied by the opposite opinions (Fig. [T]). In the case of the ring, 
an interface separates a domain of Os and a domain of Is. Figures [2](a) and Mjo) represent 
an example time course of m and En for the exponential and power-law voter models on the 
ring with N = 100, respectively. Because two interfaces that meet on a link annihilate each 
other and decrease E^ by two and because new interfaces are not produced inside a domain 
containing a single opinion, E^ monotonically decreases for both the models. Figure [2] also 
indicates that E^ = 2 for most of the time before consensus. When Ex = 2, the ring consists 
of one domain of Os and one domain of Is. Therefore, T can be approximated by the time 
at which the two interfaces collide since E^ decreases to two. For the exponential voter 
model on the ring, the distance between the two interfaces follows the simple random walk. 
Therefore, T is estimated to be the time needed for the random walker to travel a distance 



of O(iV); i.e., T oc N 2 [2j, |32j. Because the probability with which an interface moves to the 



left and that to the right are equal to 1/2 independent of the detail of p(r), the number of 




FIG. 1. Schematic of interface on the ring. The local network surrounding an interface (pointed by 
the arrow) is depicted. Open and solid circles represent voters with opinions and 1, respectively. 



(a) (b) 




time time 

FIG. 2. Time course of the fraction of 1 voters m (gray lines) and the number of interfaces En 
(black lines) in the (a) exponential and (b) power-law voter models on the ring with 100 nodes. 
The values of E^ shown are those normalized by the total number of links in the ring, which is 
equal to N. 

the movements of the interface until the consensus is reached for the power-law voter model 
is the same as that of the exponential voter model. Therefore, T for the power-law voter 
model also obeys T oc N 2 . 

Therefore, the effect of the power-law p(r) on T for the ring can be ascribed to the 
difference in the behavior of the interval between successive movements of an interface. 
We refer to this interval as the sojourn time of the interface and denote it as s. For the 
exponential and power-law voter models on the ring, T is roughly proportional to the product 
of (s) and N 2 . Therefore, the ratio of T for the power-law voter model to that of the 
exponential voter model is determined by (s). 

To derive (s) for the exponential and power-law voter models, we consider the situations 
illustrated in Figs. E2(a) and[3](b). The time lines in the figures indicate the renewal process 
on link £. We refer to the situations ) and case (b), respectively. Suppose that an 

interface moves to the right at time t. Link £ shown in Fig. [3] becomes an interface at time 
t owing to the occurrence of an event on the adjacent link. In case (a), no update event is 
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FIG. 3. Schematics of the movement of an interface on the ring. In (a), no update event occurs 
on link £ between time and time t. In (b), at least one update event occurs between time and 
time t. In (b), the time of the last update event before t is denoted as t\ ast . The nodes and links 
near link £ on the ring are depicted above the time lines of the update events. Open and solid 
circles represent voters with opinion and 1, respectively. 



assumed to have occurred on link i between time and time t. Time t is included in IEI 
Ti, which obeys pi(ri). Therefore, s = t\ — t in case (a). In case (b), at least one update 
event is assumed to have occurred on link t between time and time t. We denote the time 
at which the last update event occurs on I before time t as ti ast . Time t is included in IEI r 
that was drawn from distribution p(r) at time ti ast . Therefore, s = r — (t — ti as t) m case (b). 

We calculate (s) of the exponential and power-law voter models for cases (a) and (b). In 
case (a), the probability that T\ > r' conditioned by t is given by 



Prfr, > t' I t) - I2ll^L- ^ 
rr(T! > r | t) — frMu , )du ,- Pi(t) 



where 



Therefore, we obtain 



P 1 (u) 



Pi(u')du'. 



pi fa 1 = £ 



Pi(ti) 
Pi(t) 



Pi(t) ' 



Because T\ = s + t (see Fig. E](a)), we obtain 



Ps(s\t) 



Pi(s + t) 



For the exponential voter model, the substitution Pi(ri) = Aexp(— Ati) in Eq. 



(6) 
(7) 

(8) 

(9) 
yields 



p s {s\t) = Aexp(— As). 



(10) 



7 



Therefore, s obeys the same exponential distribution as that obeyed by r, and we obtain 



(sr p = (r), 



(11) 



where (s) denotes the mean of s with respect to the density p s (s\t). For the power-law voter 
model, the substitution pi(ri) = [(ti + c) /c]~ ( - Q ~ 1 ' ) in Eq. ([§]) yields 



Ps(s\t) 



a-2 
t + c 



s + t + c 
t + c 



>-i) 



(12) 



This conditional density is of the same form as Pi(t\) for the power-law voter model given 
by Eq. fl^D, with c in Eq. §5§ replaced by t + c and the exponent a — 1 unchanged. Therefore, 
we obtain 

+oo (a < 3), 

«r> + ^) (3 < a). 
On the basis of Eqs. (TTT1) and ffT3]) . we obtain 



power 



(13) 



(s) cxp < (s) 



power 



(14) 



in case (a), regardless of the value of a. 

In case (b), the probability density of s conditioned by t — t\ ast is given by 

/ | . . \ P(s + t- tl ast ) 

Ps{S\l — ri as tj — — " r— , 
" \}> — t last ) 

through a derivation similar to that of Eq. ([9]), where 



(15) 



P( u ) = / p(u')du. (16) 

J u 

For the exponential voter model, the substitution p(r) = Aexp(— Ar) in Eq. ffl5|) yields 



p s (s\t - ti as t) = Aexp(-As), 



(17) 



that is, s is statistically the same as in case (a) (Eq. ( II UP ) and s obeys the same exponential 
distribution as that obeyed by r. Therefore, we obtain 



(sr p = (r), 



(18) 



where (s) denotes the mean of s with respect to the density p s (s|t — ti ast ). For the power-law 
voter model, the substitution p{r) = [(a — 1) jc\ [(r + c) /c]~ a in Eq. f[T5"j) yields 



p s (s\t - ti 



ast J 



a — 1 

t — tlast + C 



S + t - tl as t + C 
^ — ^last + C 



(19) 



Equation (|T9"|) is of the same form as Eq. fll]), with c in Eq. (@]) replaced by t — ti ast + c and 
a unchanged. Therefore, we obtain 

( s }power = ^ + *-*last _ (20) 

a — 2 

Although ( s }p° wer increases with t — ti ast , it does not diverge because t — ti ast on link t is 
finite by definition. 

On the basis of Eqs. fTl8|) and (120]) . we obtain 



(s) cxp < ( s )P° wer (21) 

in case (b). 

Equations (1141) and ( 121]) predict that T for the power-law voter model is larger than 
that for the exponential voter model. The mechanism governing the enlarged consensus 
time for the power-law voter model is essentially the same as that governing the slowing 



down of epidemic dynamics in the case of long-tail IEI distributions |15l . Il6l . Il8l . 119] . T 
for the exponential voter model and power-law voter model with a = 3.5 on the ring are 
shown in Fig. IU(a). The numerical results are consistent with the theoretical prediction 

T exp < T powcr_ Wg find ypower/ycxp ^ 1 g for di ff erent va l ues Q f N ( inset Q f Fig . g]( a )). To 

understand the origin of the value of T powc 7T cxp , we note that T powcr /T cxp w (s) powc 7(s) cxp 
because most of the movements of the interfaces correspond to case (b) when t is sufficiently 
large. We approximate t — ti ast in ( s }P owcr (Eq. (1201 ) by (t — ti ast ) ~ 1.2, which we obtained 
from the direct numerical simulations of the power-law voter model on the ring. Then, the 
substitution (r) = 1 and a = 3.5 in Eqs. (HE) and flU yields T powc 7T oxp w 1.8. 

Equation (flBl indicates that T diverges for the power-law voter model with a < 3. In 
numerical simulations with a = 2.5, we confirmed that T increases with the number of runs, 
albeit slowly (Fig.BJb)). 



B. Complete graph 



In this section, we examine the exponential and power-law voter models on the complete 



graph. The average consensus time T for the exponential voter model is given by 



26) 



N 



m log — 

m 



m ) log 



1 - m 



(22) 
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FIG. 4. (Color online) (a) Average consensus time of the exponential voter model (circles) and 
power-law voter model with a = 3.5 (squares) on the ring (main panel). The solid line indicates 
T oc N 2 . For the results shown in the main panel, T poweT /T exp (squares) and (t — ti as t) f° r the 
power-law voter model (circles) are plotted against iV in the inset, (b) Average consensus time 
of the power-law voter model with a = 2.5 on the ring as a function of the number of simulation 
runs. The error bars indicate the standard deviation. We set ./V = 100. 

where mo denotes the initial value of m. Figure [5] shows the plot of T for the exponential 
and power-law voter models on the complete graph. The value of T for the power-law voter 
model with a = 2.5 and a = 3.5 are close to that for the exponential voter model. 

In contrast to the case of the ring, we cannot understand the behavior of T for the two 
voter models on the complete graph by tracking the position of a single interface. Therefore, 
we focus on the sequences of the effective events. We define the effective event as the update 
event that occurs on an inconsistent link whose two endpoints have the opposite opinions 
(essentially the same as Fig. [T]). If an update event is an effective event, magnetization m 
changes by 1/N. Otherwise, the update event does not affect m. We denote the interval 
between successive effective events on the network as r c . T is equal to the sum of r G until the 
consensus. Note that r e is generally large when there are relatively few inconsistent links. 

In Fig. El we plot the values of (M e ) and (r e ) until the consensus for the two voter models. 
We define (M e ) as the average number of the effective events until the consensus. (M e ) is 
almost the same for the exponential and power-law models and scales as iV 2 , as shown in 
Fig. E(a). (t c ) is almost the same for the two voter models and scales as N~ l , as shown in 
Fig. 0(b). Figure [6] is consistent with Fig. [5J T is equal to the product of (M e ) and (r e ). 

In the rest of this section, we examine (r c ) for the two voter models. Suppose that an 
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FIG. 5. (Color online) Average consensus time of the exponential (circles) and power-law (squares) 
voter models on the complete graph with N nodes. The solid line is the analytical solution given 
by Eq. (|22p for the exponential voter model; i.e., T = iVln2. 
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FIG. 6. (Color online) (a) Average number of the effective events {M e ) and (b) average interval 
between successive effective events (r e ) on the complete graph for the exponential (circles) and 
power-law (squares) voter models. 



effective event occurs at time t. The time to the next effective event r e is equal to the 
smallest sojourn time among those of the E\ inconsistent links. We estimate r e by the 
extremal criterion 



1 

o A 



(23) 



where Pi(r{) represents the probability density of forward recurrence time T\ of an inconsistent 

11 
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FIG. 7. (Color online) Expected interval between effective update events r e for the exponential 
and power-law voter models, i.e., Eqs. ([26]) and (J25]) . 



link and is given by 241 ] : 



Pi(Ti) = --P(r i ) 
v / 

For the exponential voter model, by substituting 



(24) 



Pi(Ti) = Aexp(-Ar ; ) 



(25) 



in Eq. fl23|) . we obtain 



_cxp 



— log [ 1 



For the power-law voter model, by substituting 



n + c 



l 



-(a-l) 



in Eq. 



we obtain 



r power = c . 



-1+1 



(26) 



(27) 



(28) 



In Fig. Tg Xp and 7-P° wer are plotted against E\. For the complete graph, E x = m(l—m)N 2 
is relatively large as compared with that for the ring. For a sufficiently large E- u both 
Eqs. (1261) and fl28|) are approximated by {t) / E\. Therefore, we obtain Tg Xp ~ r power £ or a 
large E[, which is consistent with the results shown in Fig. [7] and the numerically obtained 
(t c ) shown in Fig. E](b). The results shown in Fig. [7] are also consistent with the results for 



12 



the ring, for which En = E- x = 2 for most of the time (Fig. [2]). As shown in Fig. [3, -r power 
is considerably larger than r^ xp for a small E[. Therefore, T oc (r c )A 2 is presumably larger 
for the power-law voter model than for the exponential voter model on the ring, which is 
actually the case (Fig. HJ). 



C. Extended rings and regular random graphs 

T power /T cxp is different for the ring (Sec. IIII A[) and the complete graph (Sec. IIIIBD . 
Therefore, the effect of the power-law IEIs on x powcv /T cxp may depend on the degree of the 
node. 

To show this, we first investigate the consensus time on the extended ring with degree 
k, where k is an even number and each node is connected to up to the (fc/2)th nearest 
neighbors on each side of the conventional ring. The extended ring with k = 2 is equivalent 
to the conventional ring considered in Sec. IIII Al The extended ring approaches the complete 
graph as k — > N — 1. For the extended ring, j , p° wer y'7 lex P j s shown as a function of k by 
filled symbols in Figs. E£a) and Mb) for N = 500 and N = 1000, respectively. T powcr /T cxp 
decreases and approaches unity as k increases. 

To further examine if T power /T exp decreases with k, we investigate the consensus time 
for the regular random graph (RRG) with degree k. We generate the RRG by using the 



configuration model [33[ as follows. Each node is initially given k stubs, i.e., half links. 
Then, two stubs are chosen randomly with equal probability. We connect the two stubs to 
create a link unless a self-loop or multiple links are generated; in such a case, we discard the 
selected pair of stubs. We repeat this procedure until all the stubs are exhausted to obtain 
an instance of the RRG with degree k. If the procedure is stuck midway or the generated 
network is not connected, we restart the procedure. 

For the RRG, T powcr /T cxp is shown as a function of k by empty symbols in Figs. [8](a) and 
Mb) for A = 500 and A = 1000, respectively. 7 n P° wcr /2 lex P decreases and approaches unity 
as k increases, which is qualitatively the same as the results for the extended ring. 
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FIG. 8. (Color online) T powe 7T cxp for the extended ring (filled symbols) and for the RRG (empty 
symbols) with degree k. We set (a) N = 500 and (b) N = 1000. 



IV. CONCLUSIONS 



A. Summary of the results 



In this study, we numerically investigated voter models with the power-law IEI distribu- 
tion. On the ring, the consensus time for the power-law voter models is larger than that for 
the exponential (i.e., standard) voter model. On the complete graph, the ratio of the con- 
sensus time for the power-law voter model to that for the exponential voter model becomes 
relatively small compared to the ratio in the case of the ring. We numerically investigated 
the two voter models on the extended ring and the regular random graph, and we confirmed 
that the effect of the power-law IEI distribution to enlarge the consensus time decreases with 
the degree of the node. The difference in the consensus time originates from the fact that the 
expected time until the next update event on a link for the power-law voter model is large 
if the link has not had an update event for a long time. On the ring, the interval between 
successive movements of the interface (i.e., the link connecting nodes with the opposite opin- 
ions) is elongated by the memory effect for the power-law distribution. On the complete 
graph, the dynamics is determined by the smallest waiting time among many inconsistent 
links, which can be nearly identical for the exponential and power-law distributions. 

14 



B. Relation to previous work 



Long-tail IEI distributions are known to make epidemic dynamics slower 15|, [l6|, |l8|, Il9 |. 
Our results for the ring are consistent with these results; the occurrence of extremely large 
IEIs at some links crucially slows down the dynamics. 

Antal and colleagues considered three update rules, i.e., the link dynamics (LD), the voter 



model (VM) and the invasion process (IP) [26|, |30| (also see 



25 



31 



34J). The consensus 



time for the LD, VM, and IP is considerably different on heterogeneous networks |26l . l30 |. 
We adopted the LD update rule in this study. We can define the VM and the IP on regular 
graphs with degree k for a general IEI distribution p(r) as follows. Initially, each voter, not 
each link, is independently assigned with a random IEI until the initial update event occurs 
according to the distribution pi(ri). Suppose that an update event occurs at a voter. In the 
VM, the voter adopts the opinion of a neighbor that is selected with the equal probability 
l/k from the neighborhood. In the IP, the voter imposes its opinion on a neighbor that is 
selected with probability l/k. For either update rule, the next IEI for the voter is drawn 
from p(r). On the complete graph with N = 200, A = 1, and (a, c) = (2.5, 0.5), for example, 
we obtain T ps 1.35 x 10 2 and 3.84 x 10 4 for the exponential and power-law voter models 
with the VM update rule, respectively. With the VM, T for the power-law voter model is 
larger than that for the exponential voter model. This result is in contrast to that with 
LD; T with the LD is almost the same between the exponential and power-law voter models 
(Sec. lOB]) . We note that even on regular graphs, T with the power-law IEI distribution 
depends on the update rule, where the LD, VM, and IP are equivalent in the case of the 



exponential IEI distribution 



26 



30]. 



Finally, we remark on the relationship between our model and two other models that were 
recently proposed. Stark and colleagues examined a variant of the voter model with the VM 
update rule and an increasing inertia of voters 35j. The inertia implies that the transition 
rate at which a focal voter imitates the opinion of a neighbor decreases with the time since 
the latest change in the focal voter's opinion. They found that an appropriate amount 
of inertia shortens the consensus time on several networks. This result is opposite to our 
preliminary results for the VM described above. Although the reason of this inconsistency 
is not clear, we point out two differences in the two models. First, the transition rate can be 
infinitesimally small in the VM with the power-law IEI distribution, whereas it has a lower 



15 



bound in Stark's model. Second, when a voter experiences an update event that does not 
change its opinion, the transition rate of the voter is reset in the VM with the power-law IEI 
distribution, whereas the transition rate is not affected by such an update event in Stark's 
model. 

Fernandez- Gracia and colleagues investigated other variants of power-law voter mod- 



els 



One of their models in which the rate at which a voter experiences an update 



event decreases with the time since the last update event of the voter (called exogenous 



update in 



36]) is equivalent to the VM variant of our power-law voter model. They also 



considered the update rule in which update events occur at a constant rate (called random 



asynchronous update in 36]), which is equivalent to the VM variant of our exponential voter 
model. For the complete graph, they found that T with the exogenous update is larger than 
that with the random asynchronous update. This result is consistent with our preliminary 
result described above. 
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